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SUMMARY 


The  experimental  study  of  the  transonic  equivalence  rule  with  lift, 
previously  reported  in  Reference  has  been  extended  to  include  a  configu¬ 
ration  of  swept-back  wing-body  and  the  corresponding  equivalent  bodies  for 
lift  and  zero-lift.  For  zero  lift,  the  wing-body  is  shown  to  have  higher  wave 
drag  than  that  of  the  equivalent  body.  At  lifting  condition,  the  analyses  of 
the  data,  including  those  of  the  delta  wing-bodies  -reported  in  Reference  1, 
have  verified  the  similitude  of  drag  rise  due  to  lift  between  the  wing-bodies. 
The  additional  drag  induced  by  the  effective  area  due  to  lift  on  the  equiva¬ 
lent  body  was  found  to  agree  reasonably  with  the  wave  drag  generated  by 
lift  on  the  wing-body.  The  experiment  thus  verified  the  area  rule  with  lift 
which  must  be  considered  in  the  optimal  design  of  transonic  configurations. 


RESUME 


L ’etude  experimentale  de  la  loi  de  l’equivalence  transsonique  pour 
la  portance,  precedemment  signalee  dans  la  reference  1,  a  ete  elargie  pour  in- 
dura  le  profil  alaire  a  fleche  et  les  corps  equivalents  correspondants  pour  la 
portance  et  la  portance  nulle.  Pour  la  portance  nulle,  il  a  ete  demontre  que  le 
profil  alaire  presente  une  trainee  d’onde  superieure  a  celle  du  corps  equiva¬ 
lent.  En  condition  de  portance,  les  analyses  de  donnees,  y  compris  celles  des 
pro  fils  alaires  delta  signales  dans  la  reference  1,  ont  confirme  la  similitude  de 
l’augmentation  de  trainee  causee  par  la  portance  entre  les  divers  profils 
alaires.  La  trainee  additionnelle  induite  par  la  surface  reelle  en  raison  de  la 
portance  produite  sur  un  corps  equivalent  s’accordait  raisonnablement  a  la 
trainee  d’onde  creee  par  la  portance  sur  le  profil  alaire.  L’experience  a  ainsi 
confirme  la  loi  des  aires  pour  la  portance,  laquelle  doit  etre  prise  en  compte 
dans  la  conception  optimale  des  profils  transsoniques. 
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AN  EXPERIMENTAL  STUDY  OF  THE  TRANSONIC  EQUIVALENCE  RULE  WITH  LIFT 

PART  II 


1.0  INTRODUCTION 

In  a  previous  report  (Ref.  1)  an  experimental  investigation  of  the  transonic  area  rule  with 
lift  was  presented.  The  purpose  of  the  experiment  was  to  examine  the  validity  of  the  theoretical 
results  on  the  extension  of  the  classical  area  rule  with  lift(2).  The  extended  rule  is  given  in  a  form  of 
flow  similitude.  It  states  that  the  flows  are  equivalent  if  the  distributions  of  sources  Se  (x)  and  the 
doublets  D0  (x)  along  the  body  axis  are  identical  for  the  same  transonic  similarity  parameter  based  on 
the  maximum  cross-sectional  area  of  the  configurations.  The  sources  are  related  to  the  geometry  of  the 
body  and  its  effective  enlargement  due  to  lift,  and  the  doublets  to  the  lift  disttribution  along  the  con¬ 
figuration.  The  effective  cross-sectional  area  .’ue  to  lift  has  been  shown  to  be  appreciable  and  induces 
additional  wave  drag  in  the  transonic  flight  ranged  >. 

This  early  investigation  was  concerned  with  two  delta  wing-body  combinations  whose  drag 
characteristics  were  determined  at  transonic  speeds.  The  models  had  the  same  cross-sectional  area 
distribution  but  different  wing  span  and  thickness.  The  equivalent  bodies  for  zero  lift  and  for  a  speci¬ 
fied  lift  were  also  included  in  the  investigation.  The  delta  wing  was  chosen  because  the  planform  shape 
allowed  the  design  of  identical  effective  cross-sectional  areas  due  to  lift  for  models  with  different  swept- 
back  parameters.  However,  most  of  the  transonic  aircraft  have  swept-back  wing  with  moderate  or  large 
aspect  ratios.  From  a  practical  point  of  view,  it  is  essential  to  evaluate  the  applicability  of  the  extended 
area  rule  to  such  a  configuration.  A  swept-back  wing-body  model  was  thus  designed  to  be  investigated 
in  the  transonic  Mach  number  range  as  in  the  previous  experiment.  The  geometric  and  aerodynamic 
parameters  of  the  swept  wing  model  were  designed  to  be  closely  to  those  of  the  delta  wing,  so  that  the 
data  from  both  experiments  could  be  compared  directly. 

In  the  delta  wing  experiment,  the  wave  drag  and  induced  drag  due  to  lift  were  correlated  to 
the  square  of  lift  and  the  inverse  of  lift  curve  slope.  For  a  thin  wing  without  leading  edge  suction,  the 
induced  drag  factor,  henceforth  called  drag  factor,  defined  in  this  manner  would  approach  unity(4). 
For  the  delta  wings  with  round  leading  edges,  the  drag  factors  ranging  from  0.6  to  0.8  were  obtained. 
Based  on  the  drag  factor  defined  in  such  manner,  the  wave  drag  of  the  delta  wing  due  to  lift  was 
deduced  and  the  results  for  two  wings  were  found  to  be  closely  similar  throughout  the  Mach  number 
range  studied.  The  inclusion  of  the  lift  curve  slope  in  the  drag-lift  correlation,  however,  is  later  found 
to  be  less  desirable  as  the  former  varies  appreciable  in  the  transonic  range.  The  drag  factor  defined  in 
such  a  manner  will  no  longer  be  a  simple  correlation  function  of  the  drag-lift  polar  at  a  given  Mach 
number.  An  alternative  approach  is  to  replace  the  experimental  lift  curve  slope  by  its  theoretical 
slender  wing  value  of  7tA/2,  where  A  is  the  aspect  ratio.  The  drag  factor  is  then  in  its  conventional 
form  and  depends  only  on  the  drag-lift  correlation  at  a  given  Mach  number.  The  data  of  the  delta 
wings  in  Reference  1  and  the  present  data  of  the  swept  wing  are  analyzed  in  this  report  using  this  form 
of  the  drag  factor.  The  wave  drag  due  to  lift  is  also  deduced  in  this  form.  The  results  of  all  three  wing- 
bodies  are  then  presented  together  for  comparison. 

In  what  follows,  the  model  design  and  the  aerodynamic  consideration  of  the  swept  wing- 
body  will  be  first  presented.  The  experimental  results  and  their  analyses,  including  the  data  of  the 
delta  wing-bodies,  are  then  given  in  detail  and  related  to  the  extended  equivalence  rule  with  lift.  The 
basic  concepts  of  the  theory  can  be  found  in  Reference  1  and  will  not  be  repeated  herein. 


2.0  DESIGN  OF  THE  SWEPT  WING-BODY  MODELS 
2.1  Geometrical  Design 

The  swept  wing-body  has  a  cylindrical  central  body  similar  to  that  of  the  delta  wing-body 
Model  WB2.  The  span,  the  aspect  ratio  and  the  angle  of  sweep  of  the  leading  edge  of  the  swept  wing 
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are  taken  to  be  the  same  as  those  of  the  delta  wing  Model  WB2  before  the  wing  tip  round  off  (Fig.  1). 
The  mean  aerodynamic  chord  is  chosen  to  be  three  inches  so  that  the  wing  area  will  be  close  to  that 
of  VVB2.  The  wing  has  a  taper  ratio  of  0.4  and  the  same  symmetric  Boerstoel  profile  as  that  of  the 
delta  wings.  The  thickness  to  chord  ratio  of  the  section  is  constant  along  the  span  and  there  is  no  wing 
twist.  The  thickness  of  the  wing  is  determined  by  matching  the  cross-sectional  area  distribution  from 
the  wing  tip  forward  as  close  as  possible  to  that  of  the  delta  wing  models.  The  resulting  area  distribu¬ 
tion  for  t/'c  =  0.057  is  shown  in  Figure  2.  The  additional  area  at  the  rear  of  the  hump  comes  from  the 
outer  portion  of  the  swept  wing  downstream  from  the  leading  edge  of  the  wing  tip.  Reduction  of  the 
additional  area  to  fit  the  area  contour  of  the  delta  wing  cannot  be  achieved  without  altering  the 
geometrical  layout  of  the  planform.  It  is  preferable,  however,  to  keep  the  planform  characteristics  of 
the  swept  wing  close  to  those  of  the  WB2  wing  for  aerodynamic  consideration,  thus  no  further  modifi¬ 
cation  is  done  to  the  cross-sectional  area  distribution.  The  resulting  equivalent  body  is  thus  not  exactly 
identical  to  that  for  the  delta  wing.  The  complete  body  shape  is  shown  in  Figure  1  as  Model  BS1.  The 
wing  tip  is  rounded  off  to  ease  the  local  spanwise  loading  and  the  adverse  effects  due  to  sweep15  *.  The 
wing  tip  contour  is  designed  to  be 

v  I" x/c()  -  0.81  x/c0  >  0.8 

—  =  0.2  sin  -  (1) 

b  0.2  y/b  >  0.8 

where  b  is  the  semi-span  and  the  reference  length  c0  is  that  of  the  delta  wing  model,  taken  as  a  scale 
for  the  present  wing  design. 

The  basic  dimensions  and  the  nondimensional  parameter  of  the  final  design  for  the  wing- 
body  model  are  given  in  Table  1.  The  parameters  of  the  wing  sweep  X  and  the  thickness  rhave  values 
only  slightly  different  from  those  of  the  delta  wing  Model  WB2.  the  product  of  these  parameters  X  r 
is  nearly  identical  for  both  models.  Some  of  these  geometrical  parameters,  which  are  important  for  the 
correlations  of  the  experimental  data,  are  listed  in  Table  2  for  both  the  swept  wing  and  the  delta  wing 
models. 

2.2  Aerodynamic  Design 

When  the  wing-body  is  at  lifting  condition,  the  flow  about  the  configuration  is  governed  by 
the  effective  sources  and  doublets  distributed  along  the  body  axis.  The  strengths  of  these  singularities 
are  given  explicitly  in  terms  of  three  functions  F(x),  E(x)  and  T(x)  depending  on  the  load  distribution 
of  the  wing* 1  •-  K  For  the  present  planform,  the  load  distribution  is  unknown  and  has  to  be  determined 
experimentally  or  theoretically.  However,  it  is  noted  that  these  functions  are  not  too  sensitive  to  the 
details  of  load  distribution  for  a  given  planform.  An  approximate  expression  of  the  load  distribution 
reflecting  all  essential  characteristics  of  the  wing  loading  has  been  shown  to  be  adequate* 5  ’.  The 
expression  for  load  distributions  given  in  Reference  3  is  adopted  herein.  The  jump  of  the  velocity 
potential  across  the  wing  plane  can  be  written  as  follows, 

(x-xs)3/2 

II  0  II  =  -  aQ  +  a,  (xt  -  x)  +  a2(x,  -  x):  +  .  .  . 1 

(xt“  xt)3/2  L  J  (2) 

for  xr  <  x  <  x. 

k  i 

where  the  subscripts  t  and  t  denote  the  leading  and  the  trailing  edges  respectively.  The  wing  loading  is 
given  by  the  x-derivative  of  the  potential  jump  and  thus  gives  a  square  root  decay  of  loading  as  the 
leading  edge  is  approached.  The  coefficients  in  the  polynomial  are  determined  by  the  conditions  at  the 
trailing  edge  at  which  it  is  assumed  that  the  spanwise  loading  is  elliptical  and  the  chord-wise  loading 
decays  linearly.  In  the  wake  the  load  is  zero.  These  conditions  yield  respectively 

a„  -  -  (1  -  y2)1  /! 

1  n 


-"•*  ‘ '-r  "jt  V  .v 


With  the  conditions  at  the  trailing  edge  specified,  the  potential  jump  can  be  written  as 


(x-  x„)3/2 


0  1  =  -  (1  -  y2)1/2  -  f”  5(X  -  x)  +  2(x  -  X  )] 

n  „s/2  L  '  1  J 


(4) 


where  c  =  xt  -  xt ,  is  the  local  chord.  The  x  and  y  derivatives  of  the  potential  jump  are  given  as 
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dx.  dxn 
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5  1  dc 
2  c  dy 


3  1 


dx„ 


dy  dy 


2  x-  xe  dy  5(xt  -  x)  +  2(x  -  xc)J 


(6; 


At  the  wake  trailing  behind  the  wing,  the  conditions  in  Equation  (3)  give 


0  I  =  -  (i  -  y2)1/2 

TT 


(7 


0X  1  =  0 


(8 


*  d-y2)1/2 


(9 


matching  continuously  with  the  expressions  in  Equations  (4),  (5)  and  (6)  at  the  trailing  edge.  It  should 
be  noted  that  the  spatial  variables  x  and  y  in  this  section  (and  in  the  Appendix)  have  been  scaled  by 
the  total  length  cn  and  the  semi-span  b  of  the  wing  planform  respectively.  With  the  load  distribution 
specified,  the  functions  F(x),  E(x)  and  T(x)  are  evaluated  numerically.  Some  details  of  the  numerical 
procedure  are  given  in  the  Appendix.  The  results  are  shown  in  Figure  3  with  the  abscissa  scaled  by  the 
planform  length  cn ,  so  that  a  direct  comparison  with  the  delta  wing  case  can  be  made.  For  the  swept 
wing,  the  x-derivative  of  the  function  F(x),  Fv(x),  relating  to  the  axial  load  distribution,  peaks  at 
about  the  mid  point  of  the  planform  length,  while  on  the  delta  wing,  it  reaches  a  maximum  near  the 
trailing  edge.  The  function  T(x),  relating  to  the  load  distribution  over  the  wing  surface,  behaves 
similarly  for  both  planforms.  The  function  E(x)  signifies  the  cross-flow  kinetic  energy  and  has  the 
same  trend  for  both  cases. 


The  design  conditions  for  the  lifting  case  are  the  same  as  those  of  the  delta  wing, 
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The  effective  area  due  to  lift  can  now  be  calculated  from  the  expression  (Eq.  (25)  of  Ref.  1)  for  the 
wing  body. 
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where  the  tilded  x  is  scaled  by  the  body  length  C  and  the  barred  x  by  the  planform  length  cn.  The 
resulting  effective  area  is  shown  in  Figure  4.  The  area  distribution  for  the  delta  wing  is  also  shown  for 
comparison.  Because  of  the  nearly  symmetrical  form  of  the  loading  functions  F^(x)  and  T(x)  for  the 
swept  wing,  the  effective  area  is  more  evenly  distributed  at  both  sides  of  the  peak,  while  for  the  delta 
wing  the  distribution  is  highly  weighted  towards  the  trailing  edge.  With  the  effective  area  due  to  lift 
added  to  the  basic  geometric  cross-sectional  area  Sc(x),  the  final  area  distribution  of  the  equivalent 
body  for  the  lifting  case  is  shown  in  Figure  2.  The  corresponding  model,  designated  BS2,  is  shown  in 
Figure  1.  Note  that  the  nose  piece  is  common  for  all  models. 

3.0  EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 


3.1  Experimental  Program 

The  experiments  were  performed  in  the  transonic  test  section  of  the  NAE  5-ft  (1.5  m)  X 

5- ft  (1.5  m)  blowdown  wind  tunnel.  The  experimental  set-up  was  similar  to  that  of  the  delta  wing 
investigation.  The  forces  and  moments  of  the  model  were  measured  by  a  one  inch  Able  MK22 

6- component  internal  balance.  The  base  pressure  was  also  measured  to  complement  the  balance  drag 
measurements.  The  static  pressure  distribution  along  the  top  and  bottom  walls  were  monitored  by 
one  inch  static  pressure  pipes.  Boundary  layer  transitions  on  the  wings  and  the  nose  of  the  body  were 
left  free  as  in  the  delta  wing  experiments.  Surface  flow  visualizations  were  obtained  for  certain  Mach 
numbers  and  angles  of  attack.  Schlieren  visualizations  of  the  flow  field  were  also  taken  for  cases  with 
Mach  number  close  to  unity.  The  test  Mach  number  ranged  from  0.8  to  1.1.  The  Reynolds  number 
was  fixed  at  8  X  106  per  foot  or  2  X  106  based  on  the  mean  aerodynamic  chord  of  the  swept  wing 
model.  The  raw  data  were  reduced  by  the  same  computer  program  as  used  for  the  delta  wing  investi¬ 
gation.  The  results  are  presented  in  terms  of  C(  ,  a  and  Cp  in  Figure  5  for  the  wing-body  Model  WBS. 
The  same  data  are  also  listed  in  Table  3.  These  data  are  further  analyzed  in  terms  of  parameters  of  the 
equivalence  rule  and  compared  directly  with  those  of  the  delta  wing  models. 

3.2  Drag  at  Zero  Lift 

The  transonic  drag  rises  for  the  wing-body  Model  WBS  and  the  equivalent  body  BS1  are 
shown  in  Figure  6.  The  drag  rise  of  the  swept  wing  Model  WBS  follows  closely  to  that  of  the  delta 
wing  Model  WB2  up  to  a  value  of  K,  the  similarity  parameter  of  about  -  1.0,  corresponding  to  a 
Mach  number  of  0.98.  At  higher  Mach  numbers,  the  wave  drag  of  the  swept  wing  model  increases 
steadily  and  is  higher  than  those  of  the  delta  wing  models. 
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For  a  swept  win^  with  moderate  aspect  ratio  and  sweep-track  alible,  the  flow  about  the  wing 
retains  some  two-dimensional  sectional  characteristics  and  at  high  Mach  numbers  shocks  with  nearly 
constant  strength  extend  along  the  full  span.  This  shock  formation  can  be  seen  from  the  Schlicrcn 
flow  visualizations,  a  typical  one  shown  in  Figure  7,  taken  at  a  Mach  number  of  0.98,  shows  an  ob¬ 
lique  shock  at  the  trailing  edge  turning  rapidly  normal  to  the  flow.  For  a  delta  wing,  the  flow  about 
the  wing  is  more  three-dimensional  and  nearly  conical,  thus  the  shock  decays  along  the  span.  The 
surface  flow  visualization  in  Figure  8  shows,  at  the  outer  portion  of  the  swept  wing,  that  the  shock 
interacts  strongly  with  the  boundary  layer  causing  local  flow  separation.  The  strong  shock  and  the 
flow  separation  on  the  swept  wing  may  contribute  to  the  higher  drag  in  comparison  with  the  delta 
wing. 

The  drag  rise  of  the  equivalent  body  for  the  swept  wing-body  BSJ  is  also  shown  in  Figure  5. 
For  Mach  number  less  than  unity,  the  drag  rise  of  BS1  is  about  half  of  that  of  WBS,  although  the 
general  shape  of  these  two  curves  are  closely  similar  over  the  range  of  the  Mach  number  considered. 
The  same  discrepancy  has  been  observed  in  the  original  area  rule  study"’ 1  and  in  more  recent  experi¬ 
ments'  s  1 .  The  Schlieren  picture  of  BS1  in  Figure  7  shows  a  near  normal  shock  formed  not  far  from 
the  end  of  the  hump.  The  shock  decays  faster  than  that  of  the  WBS  and  thus  is  weaker.  The  surface 
flow  visualization  indicates  no  boundary  layer  separation  at  the  aft  portion  of  the  hump  (Fig.  8). 

Figure  (5  also  shows  that  the  drag  rise  of  BS1  is  appreciably  lower  than  that  of  Bl,  the 
equivalent  body  for  the  delta  wing  models.  Comparison  of  the  Schlieren  pictures  of  BS1  in  Figure  7 
and  that  for  Bl  in  Figure  9  of  Reference  1  clearly  indicates  that  the  gradual  reduction  of  the  hump 
at  the  downstream  side  helps  to  weaken  the  shock.  It  may  be  recalled  that  the  body  contour  ahead 
of  the  maximum  cross-section  is  practically  identical  for  both  models. 

For  the  far  field,  the  data  of  the  static  pressure  distributions  along  the  upper  and  the  lower 
walls  are  found  to  be  too  scattered  for  analysis  to  establish  flow  similarity.  Thus  no  study  of  the  far 
field  similarity  was  made  for  the  swept  wing  model  and  its  equivalent  bodies. 

:S.:i  Lift 

The  lift  versus  angle  of  attack  curves  for  all  test  Mach  numbers  are  plotted  in  Figure  5. 

Slight  nonlinearity  of  the  lift  curve  above  two  degrees  angle  of  attack  can  be  observed  for  all  cases 
especially  for  those  between  Mach  number  0,9  to  1.0.  The  nonlinear  increase  of  lift  has  been  shown  to 
be  due  to  the  tip  vortex  sheet,  which  induces  a  vortex  lift,  usually  associated  with  slender  delta 
wings"'  ’.  Another  significant  aerodynamic  chracteristic  of  the  swept  wing  observed  in  the  experiment 
is  the  strong  up  wash  over  the  outer  portion  of  the  wing,  induced  by  the  spanwise  component  of  the 
trailing  vorticity  shedded  from  the  wing' 1  0 1 .  This  effective  incrase  of  angle  of  attack  along  the  span 
causes  the  formation  of  a  strong  shock  and  flow  separation  at  the  outer  portion  of  the  wing,  which 
can  be  clearly  seen  from  the  surface  flow  visualization  ( Fig.  9).  At  2.9  degree  angle  of  attack  for 
instance,  trailing  edge  separation  starts  from  the  wing  tip  and  spreads  inward  to  about  one-third  of 
the  wing  span.  Strong  shock-boundary  layer  interaction  is  indicated  by  the  rapid  turning  of  the  surface 
streamline  at  the  separation  region.  Stronger  interaction  can  be  seen  as  the  angle  of  increase  to  5.5 
degree.  The  S  shape  of  surface  streamlines  at  the  mid  portion  of  the  wing  indicates  the  formation  of  a 
second  shock  ahead  of  the  main  shock  at  the  inner  portion  of  the  wing,  the  existence  of  which  has 
been  shown  both  experimentally  and  theoretically' 1  1  •'  2). 

The  lift  curve  slopes  at  zero  lift  are  correlated  against  the  transonic  similarity  parameter 
based  on  the  wing  thickness  in  Figure  10.  The  delta  wing  cases  are  also  shown  for  comparison.  The 
present  data  show  a  small  rise  between  Mach  numbers  of  0.92  to  0.96  with  a  peak  at  0.94.  The  rest 
of  the  data  vary  only  slightly  through  the  Mach  number  range.  The  value's  at  both  the  subsonic  and 
supersonic  ends  approach  those  of  the  delta  wing  WB1 .  At  the  low  Mach  number  end,  the  value's  e>f 
f  ’ |  .  A  agre-e  we'll  with  that  pre'die  texl  by  subsonic  line>ar  the>e>ry '  1  ' At  Mae  h  number  unity,  the-  tran¬ 
sonic  correlation  of  C|  _(t/e)i  '  versus  Aft/e')1  M  is  e-leise*  to  that  e>f  the*  delta  wing  WB2'  1  4  1 . 


3.4  Induced  Drag 


Before  the  wave  drag  becomes  appreciable,  the  drag  force  acting  on  the  wing-body  is  com¬ 
posed  of  the  induced  drag  due  to  lift,  the  skin  friction  and  the  pressure  drag  caused  by  boundary 
layer  displacement  or  separation.  The  latter  two  components  can  be  identified  as  viscous  drag  as  they 
exist  only  in  viscous  flows.  At  zero  lift,  the  total  drag  consists  of  the  viscous  drag  only.  With  the  viscous 
drag  at  zero  incidence  subtracted,  the  drag  coefficient  is  plotted  against  the  square  of  lift  coefficient 
for  Mach  number  0.8  in  Figure  11.  The  upper  and  the  lower  limits  of  the  induced  drag,  corresponding 
to  zero  and  full  leading  edge  suction  respectively,  are  also  shown  in  the  graph.  The  delta  wing  results 
from  the  previous  experiments  are  also  plotted  for  comparison.  The  drag  polar  of  the  swept  wing  is 
slightly  lower  than  that  of  the  delta  wing  WB2,  indicating  that  the  planform  shape  has  only  a  small 
effect  on  the  induced  drag  of  these  two  wings  having  nearly  the  same  aspect  ratio.  The  parabolic 
drag  polar  is  linear  up  to  C(  -  about  0.1,  above  which  the  viscous  drag  is  anticipated  and  the  correla¬ 
tion  becomes  nonlinear.  The  classical  expression  for  induced  drag  as  a  function  of  lift  and  aspect  ratio 
can  be  written  for  the  linear  portion  of  the  drag  polar  as 
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where  k  is  the  induced  drag  factor.  The  coefficients  of  the  parabolic  relation  has  been  obtained  by  a 
polynomial  fit  of  the  data  and  the  coefficient  of  the  second  degree  term  yields  k/7rA.  This  form  is 
used  in  the  present  analysis  in  preference  to  k/C,  _  employed  in  Reference  1.  The  latter  involves  the 
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lift  curve  slope  C,  _  which  varies  appreciably  in  the  transonic  range,  thus  the  drag  factor  k  is  no  longer 

a 

a  simple  coefficient  of  drag-lift  correlation. 

As  the  Mach  number  increases,  the  wave  drag  sets  in.  The  wave  drag  due  for  lift,  as  well  as 
the  induced  drag,  is  related  to  the  square  of  the  lift  forced  5  K  The  expression  in  Equation  (11)  can 
therefore  be  used  for  both  the  induced  drag,  also  called  vortex  drag,  and  the  wave  drag,  with  the 
coefficient  k  containing  both  components.  The  variation  of  k  against  Mach  number  for  the  swept  wing 
and  the  delta  wings  are  shown  in  Figure  12.  At  the  low  Mach  number  end,  the  contribution  comes 
solely  from  the  vortex  drag.  For  the  delta  wings  WB1  has  a  lower  induced  drag  factor  than  WB2, 
which  has  a  larger  aspect  ratio.  This  may  be  caused  by  the  higher  suction  generated  by  the  thicker 
wing  section  of  the  WB1  model.  As  the  Mach  number  increases,  the  drag  factor  for  the  WB1  model 
increases  to  a  maximum  and  then  drops  sharply  to  a  level  as  low  as  that  before  the  drag  rise.  For  WB2, 
the  drag  factor  peaks  at  about  the  same  Mach  number  of  0.96,  drops  slightly  and  then  increases  gra¬ 
dually  as  the  flow  becomes  supersonic.  The  swept  wing  WBS  has  the  lowest  value  of  k  before  drag  rise. 
However,  it  also  has  the  highest  value  after  drag  rise.  Before  the  drag  rise,  both  WB2  and  WBS  show  a 
dip  of  drag  factor  implying  a  reduction  of  vortex  drag.  This  slight  reduction  of  drag  is  caused  by 
favourable  interference  between  the  wing  and  the  body  of  certain  configurations  and  has  been 
observed  in  other  tests  of  wing-body  configurations* 7  ■ 1  6  \ 

3.5  Drag  Rise  of  Equivalent  Bodies  with  Lift 

The  transonic  drag  rise  of  the  equivalent  body  of  revolution  for  zero  lift  of  the  swept  wing 
configuration  BS1  has  been  discussed  in  Section  3.2.  It  was  shown  in  Figure  6  that  the  drag  rise  of  BS1 
is  much  lower  than  that  for  the  wing-body  Model  WBS.  The  drag  rise  characteristics  of  the  equivalent 
body  BS2,  with  the  effective  area  due  to  lift  forming  a  bigger  bulge,  would  be  expected  to  be  similar. 
This  indeed  is  the  case  when  both  data  are  plotted  in  terms  of  the  similarity  parameters  for  axisym- 
metrical  bodies  as  shown  in  Figure  13.  The  data  are  well  coalesced  even  though  these  two  bodies  are 
not  exactly  geometrically  similar.  The  data  of  the  equivalent  bodies  for  the  delta  wing-bodies,  B2  and 
Bl,  with  and  without  lift  respectively  are  also  shown  in  the  same  figure.  Again,  the  data  are  well 
correlated.  The  lower  drag  rise  of  the  BS1,  BS2  bodies  in  comparison  with  the  Bl,  B2  bodies  is  due  to 
broadening  of  the  hump  on  the  formers  (see  Fig.  2)  that  allows  a  more  gentle  recompression  at  the 
rear  of  the  hump.  The  resulting  rear  shock  is  comparatively  weaker  than  that  of  the  Bl,  B2  bodies 
and  no  boundary  layer  separation  has  been  observed  on  the  BS1,  BS  bodies. 


4.0  LIFT  DEPENDENT  WAVE  DRAG  OF  THE  SWEPT  WING  AND  DELTA  WING-BODIES 

The  lift  depending  wave  drag  which  is  the  centre  of  the  present  study  will  be  analyzed  in 
this  section.  The  data  for  the  delta  wing-bodies  from  Reference  1  will  also  be  considered,  together 
with  the  present  swept  wing-body  data. 

4.1  Correlations  of  the  Drag  Rise  of  the  Wing-bodies 

The  wave  drag  at  lifting  condition  can  be  deduced  by  subtracting  the  vortex  drag  portion 
from  the  total  drag.  Since  the  vortex  drag  does  not  change  with  Mach  number,  its  value  at  Mach 
number  0.8,  the  lowest  Mach  number  of  the  experiment,  is  taken  to  be  the  vortex  drag.  No  wave  drag 
is  present  at  this  Mach  number.  The  wave  drag  factor  is  therefore 

kw  =  k  -  kv  (12 

where  kv  is  the  vortex  drag  factor  and  kw  the  wave  drag  factor.  The  wave  drag  data  is  now  expressed 
in  the  same  form  as  the  vortex  drag  in  terms  of  lift  and  aspect  ratio  as 


The  resulting  wave  drag  factor  for  the  swept  wing  and  the  delta  wing-bodies  are  shown  in  Figure  14. 

It  should  be  noted  that  the  drag  coefficient  in  Equation  (13)  is  based  on  the  gross  wing  area. 

Since  the  lift  is  mainly  carried  by  the  wing,  the  scaling  length  for  the  drag  rise  due  to  lift 
should  therefore  be  the  overall  chord  c0  of  the  wing  planform(  1  7  K  The  wave  drag  rise  written  in  terms 
of  the  drag  coefficient  Cu  based  on  wing  area  is  thus 
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where  Sc  is  the  maximum  cross-sectional  area  of  the  model.  Substituting  the  wave  drag  coefficient  in 
th  .orm  of  Equation  (13)  from  the  experimental  data  into  Equation  (14),  the  drag  rise  takes  the 
following  form 


The  expression  can  be  further  written  in  terms  of  the  lift  parameter  e  —  ,  which  control  both  the 
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source  and  the  doublet  strength  due  to  lift  (see  Eq.  (10)  and  Eqs.  (12),  (13)  of  Ref.  1), 
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Since  all  wing-body  models  have  the  same  maximum  cross-sectional  area,  thus  7  A  is  the  same  for  all 
cases.  If  the  drag  rise  is  directly  related  to  the  lift  parameter,  then  the  correlation  of  the  wave  drag 
factor  as  shown  in  the  bracket  of  Equation  (16)  must  be  a  function  of  Mach  number  only.  The  corre¬ 
lations  in  such  form  for  the  data  of  three  models  as  a  function  of  Mach  number  are  shown  in  Figure  15. 
The  data  for  the  swept  wing  Model  VVBS  and  the  delta  wing  Model  WB2  collapse  well  through  the 
whole  Mach  number  range.  These  two  models  have  nearly  the  same  span  but  different  root  chords  and 
planform  shapes  which  seem  to  have  little  effect  on  the  correlation. 

The  wave  drag  factor  for  the  delta  wing  Model  VVB1  with  a  short  span  does  not  follow  the 
correlations  of  the  other  models  below  Mach  number  0.98.  The  correlation  has  much  higher  values 
than  the  other  data.  For  Mach  number  greater  than  0.98,  the  correlation  drops  to  the  same  level  as 
the  others  and  then  tends  to  increase  again  at  the  end  of  the  experimental  reach  number  range.  This 
shows  that  the  drag  rise  for  wings  with  different  spans  cannot  be  compared  for  the  same  lift  parameter 
in  which  the  lift  is  scaled  by  the  total  planform  length.  It  would  be  more  appropriate  in  the  drag  rise 
region  if  the  lift  is  caled  by  the  span.  This  form  is  readily  given  by  Equation  (15)  and  the  wave  drag 
factor  correlated  in  such  a  manner  is  shown  in  Figure  16.  The  data  of  the  delta  wing  WB1  now  co¬ 
alesce  much  better  in  the  drag  rise  region  with  the  data  of  the  other  two  wings  with  larger  spans.  This 
correlation  implies  that  the  drag  rise  due  to  lift  is  directly  related  to  Ct  -  / A  in  a  form  similar  to  that 
of  the  vortex  drag. 

From  the  Schlieren  flow  visualization  taken  during  the  experiment,  the  distinct  character¬ 
istics  of  the  flow  with  respect  to  these  two  regions  of  correlation  can  be  observed.  At  the  drag  rise 
region,  the  shock  waves  formed  about  the  wing  are  situated  ahead  of  the  trailing  edge.  The  strength  of 
the  shock  for  this  type  of  transonic  flow,  with  a  supersonic  pocket  terminated  by  a  shock,  depends 
strongly  on  the  sectional  characteristics  of  the  wing.  The  dependency  on  the  sectional  geometry  can  be 
readily  seen  from  the  fact  that  for  the  present  models  having  the  same  cross-sectional  area,  the  corre¬ 
lation  factor  Ar-  in  Equation  (15)  can  be  shown  to  be  proportional  to  the  mean  thickness  of  the 
wing.  As  the  Mach  number  increase  to  0.98  and  above,  the  shock  moves  to  the  trailing  edge  and  the 
flow  over  the  wing  is  fully  supersonic  except  at  a  small  region  around  the  leading  edge.  The  charac¬ 
teristic  of  the  wave  drag  is  then  changed  and  follows  the  classical  results  of  wave  drag  for  supersonic 
flow. 

It  may  be  recalled  that  the  wing-body  VVB1  is  designed  to  have  the  same  effective  cross- 
sectional  area  due  to  lift  as  the  Model  WB2  by  manipulating  the  special  relations  of  the  lifting  function 
F,  T  and  E  for  the  delta  planform  so  that  the  flow  similitude  may  be  established  (see  Eq.  (10)  and 
also  Section  3.1  of  Ref.  1  for  details).  The  similitude  of  flow  field  does  seem  to  exist  for  the  far  field 
as  shown  in  Figure  21  of  Reference  1.  However,  the  similitude  of  drag  rise  in  the  form  of  Equation  (16) 
does  not  seem  to  be  applicable  and  the  alternative  form  Equation  (15),  should  be  used. 

We  may  now  summarize  the  conditions  of  application  of  the  area  rule  with  lift.  For  configu¬ 
rations  with  nearly  the  same  sweep-back  parameter  A,  the  drag  rise  is  directly  related  to  the  square  of 

the  lift  parameter  c  — .  For  configurations  with  large  difference  in  A,  the  drag  rise  is  correlated  to  the 
square  of  the  lift  coefficient  and  the  aspect  ratio  as  C,  -  /A.  When  the  flow  over  the  wing  is  supersonic, 
the  lift  parameter  e  —  should  be  used  regardless  of  the  value  of  A.  The  planform  shape  has  only  minor 

T 

effect  on  the  drag  rise. 

4.2  Comparison  with  Equivalent  Bodies 

The  drag  rise  of  the  equivalent  body  at  zero  lift  has  been  shown  to  be  much  lower  than  that 
of  the  wing-body  for  the  swept  wing  experiment  (Section  3.2).  A  direct  comparison  of  the  drag  rise 
of  the  equivalent  body  including  lift  with  the  wing-body  at  lifting  condition  will  therefore  not  be 
meaningful.  However,  the  differential  of  the  drag  rise  between  these  two  equivalent  bodies  is  caused 
by  the  additional  effective  area  due  to  lift.  Thus  a  comparison  of  the  drag  rise  differential  due  to  lift 
as  given  in  the  proceeding  section  is  justified  and  will  provide  an  assessment  of  the  range  of  application 
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of  the  theory.  For  comparison,  the  drag  rise  is  now  given  in  terms  of  the  drag  coefficient  based  on  the 
maximum  cross-sectional  area,  which  is  the  same  for  the  wing-bodies  and  the  equivalent  body  without 
lift.  For  the  wing-body  the  drag  rise  is  calculated  from  the  experimental  correlations  given  in  Figures 

15  and  16  at  the  design  lift  condition  with  e  —  equal  to  0.2024. 

r 

For  the  swept  wing-body  the  comparison  is  shown  in  Figure  17.  The  drag  rise  differential 
agrees  very  well  with  the  drag  rise  to  lift  for  Mach  number  below  0.96.  For  Mach  number  greater 
than  0.96,  the  drag  rise  due  to  lift  of  the  wing-body  is  lower  than  the  drag  rise  differential  of  the 
equivalent  bodies.  A  similar  comparison  for  the  delta  wing  models  is  given  in  Figure  18.  For  the  Model 
WB1,  the  drag  rise  up  to  Mach  number  of  0.96  is  evaluated  with  CL  2  / A  equal  to  0.0078  corresponding 

to  e  — ;  of  0.0829.  For  the  Model  WB2,  the  value  of  C,  2 /A  yields  e  —  of  0.2024  the  design  value.  If 
r  r 

the  drag  rise  is  compared  with  the  same  e  —  in  this  region,  the  drag  curves  will  look  like  the  wave  drag 

T 

factor  curves  in  Figure  15.  The  comparison  in  Figure  18  shows  that  for  delta  wings  the  equivalent 
body  differential  gives  a  much  higher  drag  rise  than  that  generated  by  lift.  However,  it  should  be  noted 
that  flow  separation  occurs  at  the  rear  of  the  hump  of  the  equivalent  body  with  lift,  B2  at  the  higher 
Mach  numbers.  This  undoubtly  will  give  higher  drag  differential  in  the  ideal  case  without  separation. 
From  these  comparisons  we  may  observe  that  the  additions  drag  induced  by  the  effective  area  due  to 
lift  on  the  equivalent  body  is  close  to  the  wave  drag  generated  by  lift  in  the  drag  rise  region.  At  the 
upper  end  of  the  transonic  range  the  differential  drag  of  the  equivalent  bodies  is  higher  than  that 
generated  by  lift. 

5.0  CONCLUDING  REMARKS 

The  experimental  study  of  the  transonic  equivalence  rule  with  lift  presented  in  Reference  1 
has  been  extended  to  include  a  swept-back  wing-body  and  the  corresponding  equivalent  bodies  of 
revoluation  for  lifting  and  non-lifting  cases.  The  results  of  the  analysis  of  the  wave  drag  and  drag  rise 
in  the  transonic  range  are  summarized  as  follows: 

1.  For  zero  lift  drag  the  swept  wing-body  has  the  same  drag  rise  characteristics  as  the  delta 
wing-body  with  the  same  sweep-back  parameter  X  up  to  Mach  number  about  0.97.  At  higher 
Mach  number,  the  drag  of  the  swept  wing-body  is  higher  than  that  of  the  delta  wing-body. 

2.  The  equivalent  body  for  the  swept  wing-body  at  zero  lift  has  much  lower  wave  drag  than 
the  wing-body.  The  high  wave  drag  of  the  wing-body  is  attributed  to  the  shock  formation 
along  the  span  and  the  flow  separation  at  the  outer  portion  of  the  wing. 

3.  The  drag  rise  due  to  lift  is  directly  related  to  the  square  of  the  lift  parameter  e  —  for  con- 

-w  r 

figurations  with  the  same  sweep-back  parameter  X.  For  configurations  with  large  difference 
in  X,  the  lift  parameter  is  in  the  form  of  CL  2  /A.  When  the  flow  over  the  wing  is  supersonic, 

the  lift  parameter  e  —  should  be  used  regardless  of  the  value  of  X.  The  planform  shape  has 
r 

only  small  effect  on  the  drag  rise. 

4.  The  additional  drag  induced  by  the  effective  area  due  to  lift  on  the  equivalent  body  agrees 
reasonably  with  the  wave  drag  generated  by  lift  on  the  wing-body. 

From  these  results,  we  may  conclude  that  the  area  rule  for  zero  lift  does  not  work  well  for 
swept-back  wing-bodies  with  moderate  aspect  ratio  as  observed  by  other  authors  in  area  rule  experi¬ 
ments'6  8  K  The  area  rule  due  to  lift  as  stated  in  (3)  and  (4)  above,  works  reasonably  well.  It  is  therefore 
essental  to  include  the  additional  effective  area  due  to  lift  in  the  consideration  of  area  rule  for  the 
design  of  transonic  configurations. 
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TABLE  1 


DIMENSIONS  AND  PARAMETERS  OF 

Body  length 
Body  radius 
Semi-span 

Mean  aerodynamic  chord 
Root  chord 
Taper  chord 

Leading  edge  swept  angle 
Aspect  ratio 
Wing  thickness 
Wing  swept  parameter 
Thickness  parameter 


THE  SWEPT  WING-BODY  MODEL  WBS 

t  =  17.5  in.  (444.5  mm) 
r  =  0.875  in.  (22.23  mm) 
b  =  7.56  in.  (192.0  mm) 
cm  =  2.98  in.  (75.7  mm) 

Cr  =  4.2857  in.  (108.8  mm) 

X  =  0.4 
A  =  38.435° 

A  =  5.04 
t/c  =  0.057 
X  =  0.4320 
r  =  0.0275 


TABLE  2 

SOME  GEOMETRIC  PARAMETERS  OF  THE  SWEPT  WING  AND 
THE  DELTA  WING  MODELS 


WB1 

cn  (in.) 

6.000 

b  (in.) 

4.044 

Sw  (in.2) 

25.155 

Sc  (in.2) 

3.647 

co/c 

0.343 

A 

2.600 

X 

0.231 

r 

0.0497 

rX 

0.0115 

WB2 

WBS 

6.000 

7.715 

7.286 

7.560 

45.278 

45.106 

3.647 

3.641 

0.343 

0.441 

4.690 

5.068 

0.416 

0.432 

0.0276 

0.0275 

0.0115 

0.0119 

TABLE  3 

AERODYNAMIC  COEFFICIENTS  OF  SWEPT  WING-BODY  MODEL  WBS 


MACH  MO. 
B=/FT(  M) 

=0.799 
=  7.79 

M  =  3  1 

MACH  NO. 
0=/FT(M) 

=0.847 
=  7.76 

M=31 

AL  OH  AB 

CL  B 

CO 

AL°H  AB 

CLB 

CD 

3. 0BA4 

-0. 2378 

0.  01  71 

-2  .967* 

-0.2309 

0.0170 

3.7763 

-0.2132 

0.0157 

—  2  •  *  7  04 

-0.2079 

0.0158 

•3  .4657 

-0 .1390 

0  .0147 

-2.3439 

-0. 1 »23 

0. 01  4  7 

2.2431 

-0.1 739 

0.0139 

-2.  1  368 

-0. 1 692 

0.0140 

•1  .9400 

-0.150  2 

0  .0  131 

-1 .83  09 

-0 .1 46  1 

0.0132 

1 .6465 

-0.  1  295 

0. 01 24 

-  1  .5259 

-0.1248 

0.0126 

1  .3655 

-0.1034 

0.0119 

-1 .2446 

-0.0993 

0  .0120 

1  .0663 

-0 .0767 

0.0116 

-  C.  9469 

-0.0715 

0.01  16 

0.  7*80 

-0. 0614 

0. 01  1  3 

-  0 .66  28 

-0.0508 

0.0113 

0.5167 

-0.0441 

0.0111 

-0.3618 

-0.0231 

0.0112 

0.2  085 

-0.0130 

0.0111 

-0.  0643 

0.0050 

0. 01 1  1 

0.0750 

0.0058 

0. 01 1  1 

0.2016 

0.0153 

0.011 1 

0.3352 

0.0185 

0  .01 1 3 

0 . 4  93  4 

0.0389 

0.0112 

0.6461 

0 .0481 

0.0113 

0.  7783 

0.0625 

0.  01  13 

0.9163 

0.0681 

0.0116 

1  .0520 

0.0791 

0.01 1 7 

1 .2293 

0.1014 

0  .01 1 7 

1  .3390 

0  .0991 

0.0121 

1  .5  098 

0.1182 

0.0123 

1  .6354 

0.1268 

0.0127 

1  .7883 

0. 1 428 

0.0130 

1  .9409 

0.1591 

0.0134 

2.0789 

0 .1  647 

0.0139 

2. 2307 

0. 1 780 

0.0143 

2.3840 

0. 1956 

0. 01  49 

2  .5354 

0.2153 

0.0155 

2.6639 

0.2186 

0.0162 

2.3222 

0  .2370 

0.0168 

2 . 95  99 

0 .2449 

0 .01 76 

3.  1 099 

0. 2636 

0.01 34 

3.2312 

0. 2627 

0.  0191 

3.4094 

0.2Q1 2 

0.0202 

3.5179 

0 .2338 

0.0210 

3. 7043 

0.3158 

0 .0223 

3.  7939 

0. 3006 

0.  0228 

4  .0002 

0.3453 

0. 0246 

4.0823 

0.3227 

0  .024R 

4.2979 

0.3745 

0 .02^3 

4  .3676 

0 .3433 

0.0275 

4.5702 

0  .3893 

0.0298 

4  .6769 

0.3830 

0.0306 

4.8757 

0.4225 

0.0334 

4.9613 

0.4056 

0  .0336 

5.1695 

0.4503 

0 .037 \ 

5.2585 

0.  4275 

0. 0369 

5. 44Q9 

0. 471 5 

0.0407 

5.4968 

0.4477 

0.0  40  0 

5  .7038 

0.5029 

0  .0456 

MACH  NO1. 

95 /^T(  M» 

II  II 

.  . 

t>  o 

M=31 

MACH  NO. 
PE/FTl  M) 

=0.91 2 
=  7.  70 

M=3  1 

AL°HAB 

CL  8 

CO 

ALPH AB 

CLB 

CD 

■3.0891 

-0.271  9 

0.0186 

-3.1551 

-0.2959 

0.0203 

2.7755 

-0.2362 

0.0167 

-2.  81  19 

-0. 2554 

0.01 «2 

2  .4570 

-0.2013 

0.0152 

-2  .4990 

-0.2243 

0.0163 

■2.2271 

-0.1312 

0.0143 

-2.2503 

-0  .1986 

0.0152 

1. 9039 

-0.  1  51  7 

0.  0134 

-  1.9402 

-0.  1699 

0. 01 41 

■1  .6057 

-0.1282 

0.0128 

-1 .61 03 

-0.1378 

0.0131 

1.3313 

-0.1 1 47 

0.0122 

-1  .3301 

-0.1116 

0  .0126 

■1 . 0573 

-0,0937 

0.  01  1  8 

-  1  .05  74 

-0.0899 

0.0120 

■0.7598 

-0.0633 

0.0114 

-0 .7506 

-0.0637 

0  .01  1  7 

•0.4460 

-0  .0361 

0.0112 

-  0.4627 

-0.0411 

0. 01  1  5 

0 . 1 4  49 

-0.0075 

0.0111 

-0 .1567 

-0.01 48 

0.0114 

0.1384 

0.0115 

0 .0 1 1 3 

0.1318 

0.0054 

0.0114 

0.4  064 

0.0277 

0.01  13 

0.4357 

0.0396 

0.01 1 4 

0  .7210 

0.0590 

0.0116 

0.7206 

0.0618 

0  .01 16 

1 .0060 

0  .0309 

0.0119 

1 . 0279 

0  .0°l  3 

0.0119 

1. 3039 

0.  1  1  04 

0. 012  2 

1  .3065 

C. 1 1 21 

0.0125 

1.6029 

0.1391 

0.0128 

1 .5985 

0.1381 

0.0133 

1 .8920 

0.1534 

0.0137 

1 . 8953 

0.  1626 

0. 01 42 

2.2026 

0.1921 

0. 0149 

2 .21 04 

0.2022 

0.0154 

2.5132 

0.2311 

0.0163 

2.5230 

0  .2339 

0.0171 

2 . 8!  97 

0  .2623 

0.0181 

2. 8464 

0.  2"79  1 

0. 01  «3 

3.1044 

0.2894 

0.0199 

3.1-*56 

0.3069 

0.021 7 

3.4207 

0.3267 

0  .0227 

3.4293 

0.3353 

0  .0244 

7. 7393 

0 . 36  3  5 

0.0260 

3. 7356 

0. 3637 

0.  02 T 8 

4.0430 

0.3990 

0.0293 

4  ,03-*l 

0.4034 

0.0314 

4.3227 

0.4213 

0.0326 

4.3468 

0.4403 

0  .0  359 

4.6165 

0.  4486 

0.0366 

4  ,6533 

0.4758 

0.0411 

4.9165 

0.4769 

0  .041  1 

4 .9552 

0.5063 

0  .046! 

5.2265 

0.5160 

0  .0469 

5.2515 

0.5361 

0.0511 

5 .5 150 

0.5368 

0.0517 

5  .5468 

0.561 7 

0.0564 

5.7660 

0.5628 

0.0568 

5.8009 

0.5745 

0.0604 

tJI  Jl  Jl  P  P  P  ( J  '  jJ  t  j  u|  r\)  rv>  >\) 


TABLE  3 

AERODYNAMIC  COEFFICIENTS  OF  SWEPT  WING-BODY  MODEL  WBS  (Cont'd) 


MACH  NO. =0.031  M=3l 

pcr/CT  (M  )  =  7.68 


AL  P  H  4  9 

CL8 

CD 

-3. 0465 

-0  .2895 

0  .0226 

-  2.  7221 

-0. 2493 

0.  01  98 

-7.3914 

-0.2130 

0.0176 

-2.1  6  97 

-0  .1  952 

0.0165 

-  1  .8557 

-0.  1703 

0.0151 

-l .5415 

-0.1403 

0.0140 

-1. 2720 

-0. 1 1 75 

0.0133 

-0.9827 

-0.0975 

0.0126 

-  0.6490 

-0  .0571 

0.0120 

-0.3526 

-0.0288 

0.0119 

-0.0920 

-0.0182 

0.0118 

0.2254 

0.0203 

0.0118 

0  .5235 

0.0480 

0.0120 

0.8216 

0.0777 

0.0124 

1.1157 

0. 1050 

0.  01  28 

1 .3975 

0. 1  266 

0.0136 

1 .7203 

0.1647 

0.0148 

2.  0267 

0.1976 

0.0160 

2.3314 

0.2294 

0 .0 1 78 

2.6285 

0  .2597 

0  .0199 

2.9398 

0.2944 

0.0223 

3.2402 

0.3287 

0  .0252 

3.5350 

0.3611 

0.0235 

3.8440 

0.3973 

0.0323 

4.1334 

0.4223 

0.0361 

4.  4366 

0. 4532 

0.  04  05 

4 .7507 

0.4903 

0.0458 

5.0416 

0  .5033 

0 .0500 

5.  29  16 

0.51 86 

0.0532 

5.6262 

0.5606 

0 .0602 

5.83  84 

0  ,5751 

0.0636 

MACH  N0.= 

0.965 

£ 

II 

*-* 

Cr/FT(M)=7,73 

AL°HAB  CUR  CO 


MACH  NP. =  0.948  M=3  1 

Ro/ft(mi=7.74 


alphas 

CL  B 

CD 

-3.1560 

-0.3018 

0.0262 

-2  .8580 

-0.271 8 

0.0236 

-2 .5274 

-0.2344 

0.021 0 

-2.  2882 

-0.2052 

0.  01  94 

-1  .9651 

-0.1789 

0.0178 

-1 .6644 

-0.1 524 

0  .0165 

-  1.3P46 

-0.  1243 

0.  0154 

-1  .0989 

-0. 1065 

0.0146 

-0.7821 

-0  .0735 

0 .0139 

-  0.501 2 

-0.0524 

0.0135 

-0.1972 

-0.0219 

0.0137 

0. 0966 

0.0075 

0.01 33 

0  .3921 

0.0347 

0.0135 

0.7013 

0  .0666 

0.0138 

l.  0024 

0. 0054 

0. 01 44 

1.2867 

0.1196 

0.0152 

1 .5862 

0.1511 

0.0164 

l ,8948 

0.1826 

0.0177 

2.1783 

0.2078 

0.0193 

2 . 4  93  0 

0  .2446 

0.0216 

2.7995 

0.2750 

0.0241 

3.0913 

0.3031 

0.0267 

3.4044 

0.3397 

0.0301 
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FIG.  1:  THE  SWEPT  WING-BODY  MODEL  AND  ITS  EQUIVALENT  BODIES 


FIG.  3:  FUNCTIONS  DUE  TO  LIFT,  E(x),  T(x)  AND  Fx2  FOR  SWEPT-BACK  AND 

DELTA  PLANFORMS 


FIG.  4:  DISTRIBUTIONS  OF  EFFECTIVE  AREA  DUE  TO  LIFT  FOR  SWEPT-BACK  AND 
DELTA  PLANFORMS  AT  THE  DESIGN  LIFT  CONDITION 


EQUIVALENT  BODY  8SI,  a  =-0.05 


WING-BODY  WBS,  a  =0.01° 

FIG.  7:  SCHLIEREN  FLOW  VISUALIZATION  FOR  WING-BODY  WBS  AND  EQUIVALENT 
BODY  BS1  AT  ZERO  ANGLE  OF  ATTACK  (NOMINAL),  M^  =  0.98 


VARIATIONS  OF  LIFT  CURVE  SLOPE  AT  ZERO  LIFT  WITH  TRANSONIC  SIMILARITY  PARAMETER  BASED  ON 
THE  WING  THICKNESS  TO  CHORD  RATIO  FOR  SWEPT  WING  AND  DELTA  WING  MODELS 


SWEPT  WING-BODY  WBS 


FIG.  11:  VARIATIONS  OF  CD  WITH  CL2  FOR  SWEPT  WING  AND  DELTA  WING  MODELS  AT  MACH  NUMBER  0.8 


VARIATIONS  OF  DRAG  FACTOR  k  WITH  MACH  NUMBER  FOR  SWEPT  WING  AND  DELTA  WING  MODELS 


EQUIVALENT  BODY  OF 


EQUIVALENT  BODIES  BS1,  BS2  FOR  THE  SWEPT  WING  CASE  AND 
31,  B2  FOR  THE  DELTA  WING  CASE 


0.020 


FIG.  15:  CORRELATIONS  OF  WAVE  DRAG  FACTOR  kw  WITH  MACH  NUMBER  FOR  FIXED  LIFT  PARAMETER  e 


r' 


si 


COMPARISON  OF  DRAG-RISE  DUE  TO  LIFT  FOR  DELTA  WING-BODIES  WB1,  WB2  AND  DRAG-DIFFERENTIAL 

OF  EQUIVALENT  BODIES  B2-B1 


APPENDIX 
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EVALUATION  OF  FUNCTIONS  T(x),  E(x)  AND  F(x)  FOR  A  SWEPT  WING 


The  functions  F(x),  T(x)  and  E(x)  constitute  the  expression  for  the  effective  area  due  to 
lift  (Eq.  (10)  in  the  text).  All  involve  integrations  in  the  y-direction  and  both  T(x)  and  E(x)  involve 
double  integrals.  These  integrals  are  evaluated  numerically  as  follows. 

The  function  T(x)  has  the  form  (Eq.  (6)  of  Ref.  1) 

/  2  8  0(x,s)  Sx  11  0  (x,y)  Bx  Cn  - — r  dsdy  (Al; 


where  a]  and  a2  are  the  y-co-ordinates  of  the  leading  edges  for  each  side  of  the  wing  respectively.  The 
integrals  can  be  evaluated  in  sequency  as 

T(x)  =  —  /  2  (1  0  (x,y)  if/2  II  0x(x,s)  1  Cn  - — — 

47r  a,  Va,  'y-s 

As  x  appears  only  as  a  parameter,  the  integrands  are  now  written  without  x  and  the  square  bracket  for 
short.  For  the  integral  with  respect  to  s  the  integrand  0x(s)  is  written  for  the  integral  As  by  the  trape¬ 
zoidal  rule  as 


ds  )  dy 


(A21 


1  a  2 

T(x)  =  —  / 

4  TT 

al 


0  \  -  0X 

n  +  l  x  n 

0x(s)  =  0X  +  ~ 

x  xn  S 


-S  (S~Sn> 

n  + 1  sn  n 


0  -  0n  s  +  0n  s 

x-  n 

n  i  s 


where 


K 


0X  -  0X 

n  +  1  x  n 


S  J.1  -  S 
n+l  n 


(A3) 


and 


AS  =  Sn  +  J  -  Sn 


a2  ~  al 

N-  1 


Substituting  Equation  (A3)  into  Equation  (A2)  and  evaluating  the  integrals,  we  have 


(^xn  '  ^xs  S")  ‘l"(y)  +  0"s  !2n(y^ 


Ky)  =  - 


N  -  I 
v 


(A4 


I,„(y)  =  [(y-*n  +  ,)*n  ly-sB  +  1  I- (y- sn)  Cn  ly- sn  I  +  (Sn  +  |  -  Sn)j 


l2n(y)  =  1/2  [(y2  -  s2  +  , )  Cn  ly  -  sn+,  I-  (y2  -  s2 )  Cn  ly  -  sn  I  +y(sn+|  -  sn)  +  l/2(s 


n  +  1 


The  integral  with  respect  to  y  can  then  be  written  with  the  trapezoidal  rule  as 


T(x)  =  ^  J,  {[**(y)I(y)]m*i  +  [**(y)I(y)]  m  (j 


It  is  explicit  that  sn  #  ym  for  n  =  m. 

The  function  E(x)  has  the  form  (Eq.  (7)  of  Ref.  1) 


E(x)  =  j-  f  J  I  <t>  (x,y)  1  I  <t>  (x,s)  1  Cn  1  dsdy 
4n  aj  y  y  (y-s) 

and  can  be  evaluated  in  the  same  way  as  T(x)  with  J  <py  (x,y)  j)  replacing  ||  (x,y)  J. 

The  function  F(x)  involves  a  single  integral  (Eq.  (5)  of  Ref.  1) 


F(x)  =  /  I  0  1  dy 

al 

In  calculating  the  effective  area  due  to  lift,  the  x-derivative  of  F(x)  is  used.  Again  applying  the  trape¬ 
zoidal  rule,  the  Fx  (x)  can  be  written  as 


Fx(x)  = 


M 
2 
m  =  1 


m  +  1 


+  0v 


■i  -  yn 


+  0(x,a7 ) 


da2 

dx 


da, 


0(x, a. ) 


dx 


For  a  fixed  x  station,  the  numerical  procedure  is  carried  out  along  the  y-direction  from 
a,  (x)  to  a2  (x)  over  the  wing  and  the  wake.  The  velocity  potential  and  its  derivatives  are  given  in 
Equations  (4),  (5)  and  (6)  over  the  wing,  Equations  (7),  (8)  and  (9)  for  the  wake.  The  resulting  fune 
tions  are  shown  in  Figure  3. 
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SUMMARY/SOMMAIRE 


The  experimental  study  of  the  transonic  equivalence  rule  with  lift, 
previously  reported  in  Reference  1,  has  been  extended  to  include  a  configu¬ 
ration  of  swept-back  wing-body  and  the  corresponding  equivalent  bodies  for 
lift  and  zero-lift.  For  zero  lift,  the  wing  body  is  shown  to  have  higher  wave 
drag  than  that  of  the  equivalent  body.  At  lifting  condition,  the  analyses  of 
the  data,  including  those  of  the  delta  wing-bodies  reported  in  Reference  1, 
have  verified  the  similitude  of  drag  rise  due  to  lift  between  the  wing-bodies. 
The  additional  drag  induced  by  the  effective  area  due  to  lift  on  the  equiva¬ 
lent  body  was  found  to  agree  reasonably  with  the  wave  drag  generated  by 
lift  on  the  wing-body.  The  experiment  thus  verified  the  area  rule  with  lift 
which  must  be  considered  in  the  optimal  design  of  transonic  configurations. 
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